On the metallic conductivity of the delafossites PdCoC>2 and PtCoC>2 
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The origin of the quasi two-dimensional behavior of PdCo02 and PtCo02 is investigated by 
means of electronic structure calculations. They are performed using density functional theory in 
the generalized gradient approximation as well as the new full-potential augmented spherical wave 
method. We show that the electric conductivity is carried almost exclusively by the in-plane Pd 
(Pt) d orbitals. In contrast, the insulating C0O2 sandwich layers of octahedrally coordinated Co 
atoms may be regarded as charge carrier reservoirs. This leads to a weak electronic coupling of the 
Pd (Pt) layers. The obtained nearly cylindrical Fermi surface causes the strong anisotropy of the 
electric conductivity. 
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I. INTRODUCTION 

Transition metal oxides attract a lot of attention due 
to a great variety of physical phenomena, most of which 
go along with the ordering of some microscopic degrees 
of freedom as a function of, e.g., temperature, pressure, 
or doping. Prominent examples are the striking metal- 
insulator transitions of the vanadates^, high-T c super- 
conductivity in the cuprates, or the colossal magnetore- 
sistance observed in the manganates^ 3 -^. Cobaltates 
have aroused much interest due to the occurrence of dif- 
ferent spin states^i^. In addition, they are promising 
materials for thermoelectric application ;] 9 ' 10 . 

Known since 1873, when Friedel discovered the min- 
eral CuFeC"2, the delafossites ABO2 keep on generating 
a strong and ever increasing interes t 11 ' 12 ' 13 , especially af- 
ter Kawazoe et al. evidenced simultaneous transparency 
and p-type conductivity-!^, which laid ground for the de- 
velopment of transparent optoelectronic devices. Fur- 
thermore, the quasi two-dimensionality of the lattice and 
the triangular coordination of atoms gave rise to such 
exciting physical properties as strong anisotropies of the 
electrical conductivity and magnetic frustration effects. 

The dclafossitc structure has the space group R3m and 
results from a stacking of monoatomic triangular layers, 
see Fig.l T] 11 ' 13 . In particular, the B-atoms are at the cen- 
ters of distorted oxygen octahedra, which share edges and 
form the characteristic BO2 sandwich layers. These tri- 
layers are interlinked by linear O-A-0 bonds, resulting 
in a twofold coordination of the A-atoms. However, the 
latter have, in addition, six in-plane nearest neighbour 
A-atoms. For this reason, the structure may be likewise 
regarded as formed from single A-atom layers, which are 
intertwined by the octahedral sandwiches. For PdCoC>2 
and PtCoC>2, we will find this latter point of view partic- 
ularly useful. Finally, the oxygen atoms are tetrahcdrally 
coordinated by one A-atom and three B-atoms. Pressure 
studies reveal an increase of the structural anisotropy on 
compression indicating the high mechanical stability of 
both the octahedral sandwich layers and the O-Pd-0 




FIG. 1: (Color online) Crystal structure of PdCo02. Palla- 
dium, cobalt, and oxygen atoms are shown as green, blue, and 
red (light, big dark, and small dark) spheres, respectively. 



(O-Pt-O) dumbbells^. 

In general, interest in the delafossite-type compounds 
has concentrated quite much on the triangular arrange- 
ment of the transition-metal atoms and the resulting pos- 
sible frustration effects, which arise once localized mag- 
netic moments are established. While most of these 
oxides have been found to be antifcrromagnetic semi- 
conductors, other class members like PdCrC>2, PdCoC>2, 
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PdRh02, and PtCo02 attracted interest due to their 
rather high metallic conductivity. In particular, PdCoC>2 
has been shown to possess one of the lowest electric re- 
sistivities of normal-state oxides, even lower than that of 
Pd metal at room temperatur o n i 12 ' 16 . Yet, the conduc- 
tivity is strongly anisotropi o 11 ^ 16 . In particular, the ratio 
of the resistivities parallel and perpendicular to the c axis 
can be as large as 200 in PdCoO?^- Photocmission data 
indicate that the density of states at the Fermi energy 
can be exclusively attributed to the Pd 4d state s 12 i 13 ' 17 . 
From the combination of photoemission spectroscopy and 
inverse photoemission spectroscopy several authors con- 
cluded that the Fermi energy is located at a shallow min- 
imum of the density of states and doping might thus lead 
to rather high values of the thermoelectric powe r 16 i 17 ' 18 . 

Despite their simple chemical formulae the delafossites 
may be regarded as prototypical superlattices where the 
composition of both the A and B layers can be used to 
strongly influence the behavior of the whole system. For 
instance, in CuCrC>2, the Fermi energy falls into the Cr 
3d band, but since the Cr layers order magnetically this 
compound is a magnetic semiconductor. In contrast, as 
will be shown below, in PdCo02, the Co layers only act 
a charge reservoirs, and conduction takes place almost 
exclusively in the Pd layers. 

As a matter of fact, quite a few electronic structure cal- 
culations for delafossitc compounds have been reported in 
the literature 19 ! 20 ^ 1 ! 22 ! 23 ! 24 ^ 5 ! 26 ! 27 ! 28 . For PdCo0 2 , there 
exist linear muffin-tin orbitals calculations by Seshadri et 
al. as well as by Okabe et al£^£±. While according to the 
former authors, who also investigated PtCoC>2, the den- 
sity of states at Ep is mainly due to the Pd 4d states 
with only small contributions from the Co 3d and O 2p 
orbitals, the results by Okabe et al. are not very specific 
in this respect. For this reason, the role of these or- 
bitals for the metallic conductivity is not yet completely 
clear. In order to resolve the issue and to make closer 
connection with the photoemission data, we apply in the 
present work the new full-potential augmented spherical 
wave method to study the electronic properties of the ti- 
tle compounds. We concentrate especially on the strong 
anisotropics and on the influence of the different species 
and orbitals on the electronic properties. 



II. THEORETICAL METHOD 

The calculations are based on density-functional the- 
ory and the generalized gradient approximation (GGA)2£ 
with the local-density approximation parametrized ac- 
cording to Perdew and Wang2£. They were performed 
using the scalar-relativistic implementation of the aug- 
mented spherical wave (ASW) method (see Kefs-HHHa 
and references therein). In the ASW method, the wave 
function is expanded in atom-centered augmented spher- 
ical waves, which arc Hankcl functions and numerical so- 
lutions of Schrodinger's equation, respectively, outside 
and inside the so-called augmentation spheres. In order 



TABLE I: Experimental and calculated lattice parameters (in 
A) and atomic positions. 



compound 




a 


c 


zo 


PdCo0 2 


exp. 


2.8300 


17.743 


0.1112 




calc. 


2.8767 


17.7019 


0.1100 


PtCo0 2 


exp. 


2.8300 


17.837 


0.1140 




calc. 


2.8989 


17.458 


0.1128 



to optimize the basis set, additional augmented spherical 
waves were placed at carefully selected interstitial sites. 
The choice of these sites as well as the augmentation 
radii were automatically determined using the sphere- 
geometry optimization algorithm*^. Self-consistency was 
achieved by a highly efficient algorithm for convergence 
acceleration^. The Brillouin zone integrations were per- 
formed using the linear tetrahedron method with up to 
1469 k-points within the irreducible wedg o 33 i 36 . 

In the present work, we used a new full-potential ver- 
sion of the ASW method, which was implemented only 
very recently 3 - 7 -. In this version, the electron density and 
related quantities are given by a spherical-harmonics ex- 
pansion inside the muffin-tin spheres. In the remaining 
interstitial region, a representation in terms of atom- 
centered Hankel functions is used 3 - 8 -. However, in con- 
trast to previous related implementations, we here get 
away without needing a so-called multiple- k basis set, 
which fact allows for a very high computational speed of 
the resulting scheme. 



III. RESULTS AND DISCUSSION 

While the previous calculations were based on the crys- 
tal structure data by Prewitt et al*^, we here calculated 
these parameters from a minimization of the total en- 
ergy. To this end, in a first step for each compound the 
lattice was relaxed and after that the oxygen parame- 
ter was optimized. The results for both compounds are 
summarized in Tab. HI Note that the deviation of the cal- 
culated structural parameters from the measured ones is 
2.4% at most, which is very well within the known limits 
of the GGA and an excellent proof of the validity of the 
new full-potential ASW method. 

In order to discuss the electronic properties, we first 
concentrate on PdCoC>2 and only after that discuss the 
changes coming with the substitution of Pt for Pd. The 
electronic bands along selected high-symmetry lines of 
the first Brillouin zone of the hexagonal lattice, Fig. [21 
are displayed in Fig. [3] The corresponding partial den- 
sities of states (DOS) are shown in Fig. 0] The rather 
complicated structure of both the electronic bands and 
the DOS results from the energetical overlap of the rel- 
evant orbitals in the energy interval shown. While the 
lower part of the spectrum is dominated by O 2p states, 
the transition metal d states lead to rather sharp peaks 
in the interval from — 3.5cV to +2cV. In particular, we 
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FIG. 2: First Brillouin zone of the hexagonal lattice. 




FIG. 3: (Color online) Electronic bands of PdCo0 2 . 



recognize the ti g and e g manifolds of the Co 3c? states as 
resulting from the octahedral coordination. In represent- 
ing these partial DOS we have used a local rotated coor- 
dinate system with the Cartesian axes pointing towards 
the oxygen atoms, cr-type overlap of the O 2p states with 
the Co 3c? e g orbitals leads to the rather sharp contribu- 
tion of the latter near — 5.8eV. In contrast, due to the 
much weaker 7r-type overlap of the O 2p states with the 
tig orbitals, these states give rise to sharp peaks in the 
interval from — 2eV to Ep. Since the Fermi energy falls 
right between the tig and e g manifolds, Co turns out 
to be in a d 6 low-spin state. Our results thus provide 
further support to the picture that Co is trivalent while 
Pd is in a monovalent d 9 configuratio n 11 ' 12 ' 13 . Further- 
more, the Co and O states give only a tiny contribution 
to the electrical conductivity, which is maintained almost 
exclusively by the Pd 4c? states. The latter are further 
analyzed in Fig. [51 which displays the five Pd Ad par- 
tial DOS. Since Pd is linearly coordinated by two oxygen 
atoms parallel to the c axis and has six Pd neighbours in 
the a-b plane we used the global coordinate system to rep- 
resent these partial DOS. With this choice, contributions 
from the d xy and d x 2_ y 2 as well as from the d xz and d yz 
states arc identical. The Pd 4c? partial DOS are strongly 
influenced by the aforementioned coordination. In par- 
ticular, the six peaks of the d^ z 2_ r 2 states near —8.0, 
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FIG. 4: (Color online) Partial densities of states (DOS) of 
PdCo02- Selection of Co 3d orbitals is relative to the local 
rotated reference frame, see text. 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 





O 
Q 



Pd4d 2 2 
Pd4d 3z 2> 




-10 -8 -6 -4 -2 2 4 

(E - Bp) (eV) 

FIG. 5: (Color online) Partial Pd Ad DOS of PdCo0 2 



—6.2, —5.2, +0.6, 1.5, and +4.2 eV are complemented by 
contributions of the O 2p partial DOS of similar shape 
reflecting the strong tr-type d-p overlap along the c axis. 
In contrast, the short in-plane Pd-Pd distances of about 
2.83A (experimental value), which are only by 3% longer 
than those in metallic Pd, lead to the broad Pd d xy ^ x 2_ y 2 
bands visible in the interval from —4.5 to 2eV in Fig. [5] 
In addition, these bands give rise to the largest contribu- 
tion at Ep, whereas the d xz>vz states do not contribute 
at all. In contrast, the contribution from the d^ z 2_ r 2 or- 
bitals is rather similar to that of the d xyx 2_ y 2 orbitals. 
However, according to a k-resolved analysis, their weights 
vary on the Fermi surface. At about +0.6 eV, the con- 
tribution of the d xy x 2_ y 2 orbitals is strongly suppressed, 
and the small dispersion along the line M-K causes the 
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sharp d 3z 2_ r 2 peak seen in Fig. [5] 

The strong quasi two-dimensionality of the electronic 
states is reflected by the Fermi surface depicted in Fig. [5] 
Apart from the very small bending parallel to the c direc- 
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FIG. 7: (Color online) Partial densities of states (DOS) of 
PtCo02. Selection of Co 3d orbitals is relative to the local 
rotated reference frame, see text. 



FIG. 6: (Color online) Fermi surface of PdCo02. 

tion, it gives rise to completely in-plane Fermi velocities 
and explains the strong anisotropy in electric conductiv- 
ity. However, note that, according to the electronic bands 
shown in Fig. [3j the dispersion parallel to the direction 
T-A in general is not negligible. In passing, we mention 
the negligible dispersion along M-K, which is typical of 
tight-binding bands in a triangular lattice. These bands 
give rise to the sharp peak at about +0.6 cV in Fig. [5] 
and can thus be attributed to the Pd c? 3z 2_ r 2 states. 

Overall, our results are in good agreement with the 
previous calculations by Seshadri et al. and by Okabe et 
a L 21 i 25 . In particular, Seshadri et al. obtained a distri- 
bution of states at Ef, which is similar to ours. Further- 
more, our results are in agreement with the photoemis- 
sion data by Tanaka et al. and Higuchi et al . 12 ^ 17 , who 
attribute the metallic conductivity almost exclusively to 
the Pd 4d states. It has been argued by several au- 
thors, that the dominant contribution is due to the Pd 
d 3z 2_ r 2 orbitals, which hybridize with the Pd 5s states. 
While we find indeed a 5s contribution of the order of 
0.02 states/eV, hence, of the order of 10% of the d^ z 2_ r 2 
partial DOS at Ep, we still point to the in-plane d xy and 
d x 2_ y 2 states, which play an even greater role at Ef than 
the d% z 2_ r 2 states. As a consequence, our results demon- 
strate that the metallic conductivity is maintained by the 
in-plane d xy and d x 2_ y 2 orbitals and the in-plane part of 
the d 3z 2_ r 2 orbitals to a similar degree with a somewhat 
greater influence of the former. This finding is confirmed 
by the calculations for PtCo02, where the DOS at Ep is 
again almost exclusively due to the in-plane states. 

In general, the partial DOS obtained for PtCo02, 
which arc displayed in Fig. \7\ resemble those of the pal- 



ladium system. However, the larger extent of the Pt 5d 
orbitals as compared to the Pd Ad states leads to an in- 
creased width mainly of the in-plane d bands as well as 
an increased overlap of these states with the O 2p states. 
As a consequence, the anisotropy of the electrical conduc- 
tivity will be somewhat reduced in PtCo02 as is indeed 
observed^. 

In passing, we mention additional GGA+U calcula- 
tions for PdCo02, which however, confirmed the GGA 
results without any noticeable changes. 



IV. SUMMARY 

In summary, we have shown that the strongly 
anisotropic metallic conductivity of PdCo02 and 
PtCo02 is almost exclusively due to the Pd (Pt) layers. 
In contrast, the octahcdrally coordinated C0O2 sand- 
wiches arc insulating. As Co was found in a low-spin d 6 
configuration, these C0O2 complexes strongly suppress 
the electronic coupling between the Pd (Pt) metallic lay- 
ers, and the Co 3+ layers merely act as charge reservoirs. 
In addition, they pronounce the quasi two-dimensionality 
of the system. It may be speculated how the introduction 
of impurities into the cobalt layers, i.e. replacing, e.g. Fe 
or Ni for Co, alters the electronic properties locally and 
might thus be used to design nano-structured materials. 
Work along this line is in progress. In this context one 
might ask to what extent the structural changes induced 
by the introduction of impurities with different covalent 
radii as compared to Co might be damped by the inter- 
mediate oxygen layers. This would provide a situation, 
where chemical substitution only acts on the electronic 
degrees of freedom. 
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